Competition between attractive and repulsive interactions in two-component 
Bose-Einstein condensates trapped in an optical lattice 
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We consider effects of inter-species attraction on two-component gap solitons (GSs) in the binary 
BEC with intra-species repulsion, trapped in the one-dimensional optical lattice (OL). Systematic 
simulations of the coupled Gross-Pitaevskii equations (GPEs) corroborate an assumption that, be- 
cause the effective mass of GSs is negative, the inter-species attraction may split the two-component 
soliton. Two critical values, ki and K2, of the OL strength (k) are identified. Two-species GSs with 
fully overlapping wave functions are stable in strong lattices (k > Ki). In an intermediate region, 
k\ > k > K2, the soliton splits into a double- humped state with separated components. Finally, 
in weak lattices (k < K2), the splitting generates a pair of freely moving single-species GSs. We 
present and explain the dependence of ki and K2 on thenumber of atoms (total norm), and on the 
relative strength of the competing inter-species attraction and intra-species repulsion. The splitting 
of asymmetric solitons, with unequal norms of the two species, is briefly considered too. It is found 
and explained that the splitting threshold grows with the increase of the asymmetry. 
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I. INTRODUCTION AND THE MODEL 

Self-supporting localized patterns in Bose-Einstein 
condensates (BECs), which are frequently called matter- 
wave solitons, have been the subject of many theoretical 
and experimental works. The solitons were first experi- 
mentally created in the condensate of 7 Li atoms loaded in 
a strongly elongated ( "cigar-shaped" ) trap [H . The use of 
the Feshbach-resonance (FR) technique [2| makes it pos- 
sible to tune the scattering length of the inter-atomic col- 
lisions in the condensate to a small negative value, thus 
providing for the weak self-attractive nonlinearity, which 
is necessary for the creation of stable solitons. In another 
experiment, the solitons were observed in the 85 Rb BEC 
remaining in the trap after the onset of collapse induced 
by the switch of the interaction between atoms from re- 
pulsive to attractive Q. Potentially, the FR technique 
in combination with a quasi-one-dimensional optical lat- 
tice (OL, i.e., a spatially-periodic potential induced by a 
pair of counter-propagating laser beams) can be used to 
create solitons with a fully three-dimensional shape 

In most experiments, BEC is created in ultracold gases 
with repulsion between atoms. In this case, bright soli- 
tons were predicted as a result of the balance between the 
repulsive nonlinearity and negative effective mass of the 
matter waves, induced by the OL The correspond- 
ing gap solitons (GSs) emerge in finite bandgaps of the 
OL's linear spectrum, where the negative effective mass 
is available. Theoretical models of matter-wave GSs were 
reviewed in Ref. and their stability was analyzed in 
detail in Refs. Q Additionally, in Ref. Q, the stability 
of the GS was also studied against quantum fluctuations, 
following the lines of the approach which was earlier de- 
veloped for ordinary matter- wave solitons in Ref. [j| . Ex- 



perimentally, GSs were created in the 87 Rb condensate 
trapped in a cigar-shaped potential combined with an 
OL applied in the axial direction |10j, with ~ 250 atoms 
in the established soliton. An essential ingredient of the 
experiment was the acceleration of the condensate, with 
the aim to push the atomic waves into the spectral region 
featuring the negative effective mass. Another approach 
to the creation of the GS was proposed in Ref. one 
may add a strong parabolic trap to the OL potential, 
confining the entire condensate to a small spatial region, 
and then gradually relax the extra trap, which may al- 
low the atomic cloud to remain in a relatively compact 
GS state. Besides that, one may expect that a chain 
of GSs may develop from the modulational instability of 
nonlinear quasi-periodic Bloch waves trapped in the OL 

1- 

Binary mixtures of BECs are also available as media 
in which various matter-wave patterns may be created. 
Most typically, the mixtures are formed by two different 
hyperfine states of the same atom, such as 87 Rb [13] and 
23 Na [l3|. Creation of a heteronuclear mixture of 41 K 
and 87 Rb was also reported • As concerns the effective 
nonlinearity in the mixture, it is known that the sign and 
magnitude of the scattering lengths which characterize 
collisions between atoms belonging to different species 
may also be controlled by means of the FR technique 

m 

In view of the latter possibility, it is reasonable to con- 
sider a binary BEC with the (natural) intra-species self- 
repulsion, while the inter-species interaction is switched 
to attraction. In recent theoretical works [l6j], it was 
proposed to use this setting to create symbiotic bright 
solitons: while self-repulsive species cannot support iso- 
lated solitons, the inter-species attraction may help to 
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create two-component solitons. Moreover, a similar per- 
spective was discussed in the context of Bose- Fermi mix- 
tures, where the interaction between bosons is repulsive, 
but the bosons and fermions attract each other [13] • Sim- 
ilarly, solitons in a binary degenerate Fermi gas, sup- 
ported by the attraction between the fermion species, 
were predicted [Hj]. It was also proposed to use the at- 
traction between fermions and bosons for making bosonic 
quantum dots that can trap fermion atoms (in particu- 
lar, gap solitons in the BEC trapped in the OL may play 
the role of such dots) . 

As mentioned above, a salient property of GSs is the 
negative effective mass, which should make their dy- 
namical behavior drastically different from that of or- 
dinary quasi-particles. In particular, both one- [21] 
and two-dimensional [22j GSs are expelled by the usual 
parabolic trapping potential, while being retained by 
the anti-trapping (inverted) potential. The negative ef- 
fective mass may essentially affect the stability of two- 
component GSs. If the BEC species repel each other, 
the repulsive force, acting in combination with the neg- 
ative effective mass, may actually keep the two compo- 
nents together. In Ref. [20], this possibility was verified, 
for two- and one-dimensional GSs in the model of the 
binary BEC. It was demonstrated that, even with zero 
intra-species interaction, the repulsion between the com- 
ponents was sufficient to generate a family of symbiotic 
gap solitons. Adding nonzero intra-species repulsion ex- 
pands the stability region of such GSs. Note that the 
symbiotic GSs found in Ref. [20] could be both of intra- 
gap and inter-gap types, i.e., with the chemical potentials 
of the two components belonging to the same or different 
finite bandgaps in the linear spectrum induced by the 
OL. 

The objective of the present work is to consider effects 
of the attraction between BEC species on two-component 
GSs trapped in an OL, in the case when the intra-species 
interactions are repulsive. As mentioned above, the nec- 
essary signs and magnitudes of the respective nonlinear 
coefficients may be adjusted by means of the FR tech- 
nique. Because of the negative sign of the effective GS 
mass, one may expect that the inter-species attraction 
destabilizes the two-component GSs, with a trend to split 
them into separate single-species solitons. This expecta- 
tion is confirmed below, by means of systematic simula- 
tions of the coupled Gross-Pitaevskii equations (GPEs) 
for the macroscopic wave functions of the two compo- 
nents, ip and <p. 

In the normalized form (we set atomic mass, Planck's 
constant, and the overall nonlincarity coefficient equal to 
1, and the OL period equal to 7r), the coupled GPEs take 
the well-known form [231 ]. 

%il>t + (l/2)tfe* - [(cos0) + (sin0) |0| 2 ] V 

+k cos(2a;)'!/> = 0, (1) 

uk + (l/2)cp xx - [(cos0) |0| 2 + (sin0) H 2 ] 

+kcos(2x)0 = 0, (2) 



where k is the strength of the OL (actually measured 
in units of the corresponding recoil energy), and angle 
8 is a parameter that determines the relative strength 
and sign of the inter- and intra-species interactions. The 
case of the intra-species repulsion and inter-species at- 
traction, which we aim to consider in this work, corre- 
sponds to — 7r/2 < 8 < 0. If ip and <j> represent two 
components of a spinor BEC (therefore, equal atomic 
masses are assumed in the two equations) with oppo- 
site z-components of the hyperfine spin, mp = ±1, then 
tan# = (ao + ^2) / (^0 — ^2), with coefficients eto and 
ao accounting for the mean-field (spin-independent) and 
spin-exchange interactions between the atoms [24j . These 
coefficient may be, if necessary, controlled by means of 
the FR technique, as mentioned above. 

Equations (JTJ and @ conserve the correspond- 
ing Hamiltonian, and two norms (scaled num- 
bers of atoms in the different species), $ = 

\^{x)\ 2 , 10(a;)| 2 } dx. We will chiefly consider sym- 
metric GS complexes, with = A$ = N (Section II); 
asymmetric states with unequal norms will also be con- 
sidered, but briefly, in Section III. 

II. RESULTS: SYMMETRIC GAP SOLITONS 

Stationary soliton solutions to Eqs. ([T|) and ^ are 
looked for as 

ifj(x, t), <f>(x, t) = e - *"** (a), e-""^), (3) 

where /i and v are chemical potentials of the two compo- 
nents. In this work, we focus on the two-component GSs 
of the most fundamental type, with both fi and v falling 
in the first finite bandgap of the linear spectrum. 

The starting point is the symmetric soliton with iden- 
tical components, fi = v and ^(x) = where the 
real wave function, 4' (a;), obeys the ordinary stationary 
equation for the single-component BEC, 

+ (1/2)*" - V2 (sin (6> + tt/4)) 4» 3 + KCOs(2a;)* = 0. 

(4) 

This equation has GS solutions provided that 
sin (8 + 7r/4) > 0, which actually implies — 7r/4 < 8 < 0, 
since we are interested in 8 < (the intcr-species 
attraction), as said above. The remaining interval, 
— 7r/2 < 8 < — 7r/4, corresponds to the ordinary symbi- 
otic solitons [l6| trapped in the OL, with \x falling in the 
semi-infinite gap. 

Equation can be solved by means of known numer- 
ical methods [a, [2l| (the variational approximation [2(| , 
which was first used in the framework of the GPE in 
Ref. [25[, becomes quite cumbersome if applied to GSs 
[27l ] , therefore we do not resort to this method here) . We 
tested the stability of symmetric solitons, generated by 
Eq. (HI), against small random perturbations (which in- 
clude a small initial separation between the components) 
by means of direct simulations of Eqs. ([T]) and In 
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FIG. 1: (Color online) Three generic scenarios of the evolu- 
tion of two-component symmetric gap solitons are shown by 
means of density-level contours of the two species in the (a;, t) 
plane. a,b) Splitting into freely moving solitons in a weak 
lattice, K = 0.2. c,d) Formation of a stable bound states of 
separated solitons in a moderately strong lattice, k = 0.4. 
e,f) Stability of the bound state in a strong lattice, k = 0.6. 
Other parameters are fixed: 9 — —0.3, norm of each compo- 
nent N = 1.5, and the evolution time, t — 500. 



agreement with the expectation that the effective nega- 
tive mass of the GS can make the two-component bound 
state unstable against the splitting under the action of 
the inter-species attraction, we have observed three dif- 
ferent scenarios of the perturbation-induced dynamics, 
depending on the OL strength, k. As shown in Figs. Q] 
a) and c), the bound state splits into two freely moving 
single-species solitons in the weak lattice, and remains 
stable in the strong OL. In the intermediate case, Fig. 
Q]b), the bound states also splits, but the components 
cannot move freely; instead, they get pinned at a finite 
distance between them. As concerns the free motion of 
the single-component soliton, in Ref. 21>\ it was demon- 
strated that the GS moves without any tangible loss if 
its amplitude does not exceed a certain maximum value, 
above which the moving soliton is being braked by the 
underlying lattice. 

A typical example of the stable symmetric soliton with 
separated centers of its two components is displayed in 
Fig. In addition, the separation between centers of the 
two species in the stable solitons of this type, generated 
by the splitting of the unstable soliton with overlapping 
wave functions, 'F = is shown in[3l According to Fig. 
[TJ the separation is zero at k > 0.46, as the overlapping 
bound state is stable in that case, while at k < 0.29 
the split components do not come to a halt (i.e., the 




FIG. 2: (Color online) A stable bound state of two gap soli- 
tons formed in a moderately strong lattice as a result of the 
splitting of the unstable state with overlapping wave func- 
tions, cf. Fig. Q] b) (parameters are as in that figure, except 
that k = 0.3). 




FIG. 3: (Color online) The separation between the peaks of 
the \E r - and ^-components in the stable two-component soli- 
ton, generated by splitting of the unstable state with overlap- 
ping wave functions [see Fig. [T]b)], versus the lattice strength. 



separation is infinite). 

As suggested by Fig. [3l two critical values, n\ and k.2, 
of OL strength k may be identified, for given 9 and norm 
N: at k < Ki, the original soliton, with &(x) = 
is unstable and splits, and at n < K2 (with K2 < Hi) 
the splinters (single-component GSs) are not pinned by 
the lattice, but rather move freely. The dependence of 
both critical values on the norm of the initial overlapping 
bound state is shown in Fig. Q] (to interpret the results 
in physical terms, recall that k shown in these plots is 
measured in units of the recoil energy). The decrease 
of Ki and K2 with increase of N is easy to understand, 
because the solitons with larger N are narrower, hence 
they are stronger pinned by the lattice (2lT |. 

In addition, Fig. [5] shows the dependence of K\ and K2 
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FIG. 4: (Color online) Critical values of the lattice depth, 
Ki (solid line), below which the solitons with the overlapping 
wave functions, \& = $, start to split, and Ki (dashed line), 
below which the single-component solitons generated by the 
splitting move freely, are shown as functions of norm N of 
each component. 




FIG. 5: (Color online) The critical lattice strengths, Ki (solid) 
and K2 (dashed), as functions of angle 6, which determines 
the relative strength of the inter-species attraction and intra- 
species repulsion. 



on 9, at fixed N. The increase of the critical values of the 
lattice strength with the growth of \6\ is quite natural: 
as the inter-species attraction, which actually pushes the 
two GS components apart, gets stronger with respect to 
the intra-species repulsion, which tends to stabilize the 
GS, it becomes more difficult for the lattice to hold the 
components together. Both critical values diverge in the 
limit of \9\ — > 7r/4, since point 9 = —ir/4 is a singular 
one, as follows from Eq. ((4]). For values of 9 close to zero, 
the threshold K2 is almost constant, because it is mainly 
determined by the ability of the lattice to stop the motion 
of the single-component soliton after the splitting. 




FIG. 6: (Color online) A stable unsplit soliton, with unequal 
norms in the two components: iV* = 1.5, iV$ = 1.6, and 
9 = -0.3, k = 1. 



FIG. 7: (Color online) Spontaneous splitting of an asymmet- 
ric soliton below the threshold, ki, for AT* = 1.5, N<s> — 1.6, 
9 — —0.3, and k = 0.5 (in this case, ki is slightly larger than 
0.5, see Fig. |8j. Total evolution time is t = 1000. 



III. ASYMMETRIC GAP SOLITONS 

We have also investigated the dynamics of asymmet- 
ric solitons, with different norms of the two components, 
^ Nq>. An example of a stable soliton of that type, 
with coinciding centers of its components, is displayed in 
Fig. 1 

Similarly to the symmetric case, the asymmetric soli- 
tons become unstable against the splitting of the two 
components if the lattice strength falls below the thresh- 
old value, Ki, see an example in Fig. [7] In the simu- 
lations, the splitting follows onset of oscillations of the 
solitons, with a significant amplitude. 

Naturally, K\ depends on the asymmetry parameter, 
(JV* - A/$) I (AT* + N$), as shown in Fig. |5J The growth 
of the splitting threshold with the increase of the asym- 
metry can be readily explained. Indeed, as the norm of 
component $ becomes smaller, its width increases. On 
the other hand, the effective pinning force (the ampli- 
tude of the effective Peierls-Nabarro potential (28|) act- 
ing on the soliton exponentially decays with the increase 
of the soliton's width, hence a higher strength of the OL 
is required to prevent the splitting instability of the two- 
component asymmetric soliton. 
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FIG. 8: (Color online) The dependence of the lattice-strength 
splitting threshold, Ki, for asymmetric solitons on the norm 
difference between the components (i.e., the asymmetry de- 
gree). 

IV. CONCLUSIONS 

We have considered effects of the attraction between 
two species in the binary BEC with intra-species repul- 
sion, trapped in the OL (optical lattice), on the two- 
component GSs (gap solitons), supported by the balance 
between the repulsion and negative effective mass in- 
duced by the OL potential. Systematic simulations con- 
firm the prediction suggested by the fact that the effective 
mass of the gap soliton is negative: if the OL is not strong 
enough, the inter-species attraction results in splitting of 
the two-component GS. We have identified two thresh- 
old values, Ki and k-z, of the OL strength (k), for the 
two-component GS with equal norms of its components. 



The unsplit solitons (with fully overlapping wave func- 
tions, \& = $) are stable in strong lattices, with k > K\) 
they split into a stationary symmetric state with sepa- 
rated components in interval k\ > k > K2 , and, in weak 
lattices, with k < k,2, the splitting generates a pair of 
freely traveling single-species GSs. The dependences of 
K\ and K2 on the norm of the original overlapping soliton, 
and on the relative strength of the inter-species attraction 
and intra-species repulsion (9) were found and explained. 
We have also considered, in a brief form, the dynamics 
of asymmetric solitons, with unequal norms of the two 
species. In particular, it was found (and explained) that 
the splitting threshold, k\, grows with the increase of the 
relative asymmetry. 

The model introduced in this work calls for further 
analysis. In particular, it may be interesting to ex- 
plore effects of intra-species attraction on symbiotic GSs 
supported by the inter-species rejmlsion (such as two- 
component GSs found in Ref. [20|), i.e., the case ex- 
actly opposite to that considered above. Another rel- 
evant generalization may be to study similar effects in 
two-dimensional GSs build of two species. 
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